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Abstract 

The radiative decay B (pK^ is observed for the first time. The branching fraction for the 
charged B~ (f)K^j decay mode is measured to be B{B^ (pK^^) = (3.4 ± 0.9 it 0.4) x 10^^. 
The photon energy distribution for the B~ — > (j)K~'y decay is presented. The signal for 
the neutral B^ — > (j)K^^ decay mode is not statistically significant and an upper limit, 
B{B^ ^K^-l) < 8.3 X 10-6 at 90% CL, is set. The analysis is based on a dataset of 90 fb"^ 
collected by the Belle experiment at the e^e~ asymmetric collider KEKB. 

PACS numbers: 13.25.Hw, 14.40.Nd 
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Radiative penguin B meson decays provide an important tool to search for physics be- 
yond the Standard ModeL Recent experimental studies [H, 0, 0, 0, IE 1^ of inclusive 
and exclusive radiative B decays are in good agreement with Standard Model predictions 
j^, New information on radiative B decays is important to further test theoretical 

models. In this analysis the exclusive decay mode B — > (pRj is studied for the first time. 
The experimental techniques used in this analysis are similar to those used in the recently 
ublished studies of B^ — K~^'ir~'j and B~^ —>■ K~^7[~TT^'~f decays by the Belle collaboration 
Measurement of the branching fraction for the B (pK'j decay and its contribution to 
inclusive radiative B —>■ Xsj decays is important to constrain theoretical models. Moreover, 
due to the narrow width of the resonance, exclusive B (pK'-f decays are well-separated 
from background and can be effectively used for measurements of the photon momentum 
over a wide interval. Such measurements can shed some light on the behavior of the photon 
momentum spectrum in inclusive decays, where the theoretically interesting region below 
2 GeV/c is difficult to study experimentally because of large backgrounds. The decay chan- 
nel B^ (pK^'-f can also be used in future high statistics measurements of time-dependent 
CP violation parameters. With a larger dataset, such three-body hadronic final state de- 
cays could also be used for angular distribution measurements The decay B~ — > (j)K~^ 
(charge conjugate modes are implied throughout this letter) can be described by conven- 
tional radiative penguin diagrams with the creation of an additional ss pair (for example 
Fig. 1). 




FIG. 1: A penguin diagram for B — > (j)K 7 with ss pair creation. 
The 90 fb~^ data sample containing (95.8 ± 0.5) x 10^ produced BB pairs was collected 



with the Belle detector Ij] at KEKB jl^ , an asymmetric energy double-storage- ring collider 
with 8GeV electrons and 3.5 GeV positrons. Belle is a general-purpose large- solid- angle 
detector that consists of a three- layer Silicon Vertex Detector (SVD), a 50-layer Central 
Drift Chamber (CDC), an array of Aerogel Cerenkov Counters (ACC), a Time-of-Flight 
Counter system (TOF), and a CsI(Tl) Electromagnetic Calorimeter (ECL) located inside 
a superconducting solenoid coil with a 1.5 T magnetic field. An iron fiux- return located 
outside the coil is instrumented to identify i^^ and muons (KLM). The detector is described 
in detail elsewhere A detailed GEANT-based simulation of the Belle detector is used 
to produce Monte Carlo (MC) event samples and determine efficiencies. 

Charged tracks with impact parameters less than 2 cm radially and less than 5 cm in z 
(the z axis is antiparallel to the positron beam direction) are used. Kaon and pion mass 
hypotheses are assigned to charged tracks with a momentum larger than 100 MeV/c, using 
a likelihood ratio Ck/-ki obtained by combining information from the CDC (dE/dx), ACC 
and TOF systems. The likelihood ratio Cr/tt ranges between and 1 and is required to 
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be larger than 0.8 for the kaon candidates and less than 0.8 for the pion candidates. A 
relaxed likelihood ratio requirement Cki-k > 0.4 is applied for the kaon candidates used to 
reconstruct mesons. 

A candidate primary photon (7) is required to have an energy E* in the T(4S) center- 
of-mass (CM) frame between 2.0 and 2.7 GeV and to lie within the acceptance of the 
barrel ECL (33° < 6.y < 128°). The main background sources of high energy photons are 
77 and rj '-ff decays. To reduce these backgrounds, restrictions are imposed on 
the invariant mass of the candidate photon and any other photon (7') in the event. The 
candidate photon is rejected if 120MeV/c2 < M(77') < 145MeV/c2 and > 30MeV 
or if 510MeV/c2 < M(77') < STOMeV/c^ and > 200 MeV. To reduce the background 
from 7T^ decays when the two daughter photons form a single cluster in the calorimeter, the 
ratio of the energy deposition in 3x3 ECL cells compared to that in 5x5 cells around the 
maximum energy ECL cell is required to exceed 95%. 

The Kg candidates are formed from tt'^tt~ combinations that have an invariant mass 
within ±10MeV/c^ of the nominal Kg mass (~ 3a in the Kg mass resolution). The two 
pions are required to have a common vertex that is displaced from the interaction point 
by more than 0.5 cm in the plane perpendicular to the beam direction. The difference 
in z coordinates for the tracks constituting the secondary vertex must be less than 2 cm. 
The angle a between the Kg flight direction and the measured Kg momentum direction is 
required to satisfy cosa > 0.8. Opposite-sign K mesons are combined to form candidates; 
their invariant mass is required to be within ±10MeV/c^ (~ 3a) of the nominal value for 
the (J) mass. 

The (f)K~'-f and (pKg'-f combinations are selected to form B~ and candidates. Two 
kinematic variables are used to extract the B meson signal: the energy difference AE = 
E"^ — -Ebeam ^^^1 the beam-coustraiued mass Mbc = \/(-£'beam)^ ~ (Pb)^' where E*^ and p*^ 
are the CM energy and momentum of the B candidate and -Ebeam beam energy. 

The events that satisfy loose requirements > 5.2GeV/c^ and \AE\ < 0.4 GeV are 
selected for further analysis. 

With these selection criteria, the primary background source is continuum e~^e~ — *• qq 
production, where q may be a u, d, s, or c quark. To separate spherical BB events from 
jet-like continuum events, a Fisher discriminant is formed from six modified Fox- Wolfram 
moments 0, Q| • Signal and background probability density functions (PDF) for the Fisher 
discriminant and the cosine of the B flight direction with respect to the z axis (cos 6'^) are ob- 
tained from signal MC and sideband data. The signal (background) PDFs are multiplied to 
form a signal (background) likelihood Cg {Cbg)- The likelihood ratio CTZ = Cg/{Cg + Cbg) 
is required to be greater than 0.3. This event topology requirement retains 92% of the signal 
events while removing 55% of the continuum events. Finally, the ratio of the second to the 
zeroth Fox- Wolfram moment, calculated using all particles in the event, is required to be 
less than 0.5. 

Signal MC studies of the AE distribution indicate that the width is dominated by pho- 
ton energy smearing, and the shape is expected to be asymmetric due to photon energy 
leakage. The resolution in Mbc of 3.2MeV/c^ is dominated by the beam energy spread 
of KEKB and is slightly improved to 3.0MeV/c^ by rescaling the photon energy so that 
AE = E^+ E*j^ + E* — -Ebeam = 0' where E'^ and E]^ are meson and K meson energies 
in the CM frame. 

The Mbc distribution for the interval -0.08 GeV < AE < 0.05 GeV and AE distribution 
for the interval 5.27 GeV/ c^ < Mbc < 5.29 GeV/ c^ for B~ (f)K~'y decay mode are shown in 
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FIG. 2: The (a) Mbc, (b) AE, (c) M{K+K~), (d) and (e) M{<j)K-) distributions for the 
B" (j)K~'y decay mode and the (f) Mbc distribution for the (pK^^ decay mode are shown 

by histograms. The sohd histograms are obtained with events from the B mass sideband. The 
curves show the results of the fit described in the text. The measured M((pK^) distribution is 
compared to those obtained from the MC simulation, with a three-body phase-space model (circles 
with error bars) or adjusted to follow the data (squares connected by a line). 



Figs. 2(a) and 2(b). To avoid systematic uncertainties that could arise from the description 
of the AE distribution, the signal yield is extracted from a fit of the Mbc distribution after 
applying the asymmetric cut —0.08 GeV < AE < 0.05 GeV. The Mbc distribution is fitted 
to the sum of a Gaussian function and the so-called ARGUS background function [l^ • The 
width of the signal Gaussian is determined from MC, while the peak position is fixed to 
5.279 GeV/c^. The background shape was studied using Mbc and (j) mass sidebands and 
found to be flat. The fitted number of B candidates is = 21.6 ± 5.6. As a crosscheck, 
the AE distribution was fitted with the Crystal Ball line shape function (the shape is 
fixed from MC) to describe the signal and a linear function to describe the background. The 
signal yield is A^ = 23.6 ±5.4 events, which is consistent with the Mbc fit result. The K^K~ 
mass distribution with two combinations per event is shown in Fig. 2(c) for events in the B 
signal region, 5.27GeV/c2 < Mbc < 5.29GeV/c2 and -0.08 GeV < AE < 0.05 GeV. The 
solid histogram in Fig. 2(c) shows events from the B meson mass sideband, 5.2GeV/c^ < 
Mbc < 5.26GeV/c^, with a normalization obtained from the Mbc distribution fit. 

The photon energy distribution for the B~ (j)K~'y decay channel for events within 
the B signal region with an extended photon energy interval 1.0 GeV < E* < 3.0 GeV is 
shown in Fig. 2(d). The solid histogram shows events from the B mass sideband, which are 
used to describe the background. The MC determined B meson reconstruction efficiency is 
nearly constant for photon energies between 1.5 GeV and 2.7 GeV. Figure 2(d) demonstrates 
that the signal events are concentrated in the range 2.1 GeV < E* < 2.5 GeV, and events 
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outside that range are consistent with background, estimated from the B mass sideband 
events. The bulk of inclusive B Xg'y events were also observed in the E* > 2.0 GeV 
range 0, 0, • The maximum photon energy is kinematically constrained by the mass of 
the hadronic system and is therefore slightly lower in B~ — > (pK'j decays compared to 
B ^ X,7. 

To search for a possible contribution from kaonic resonances decaying to (f)K~, the (j)K~ 
invariant mass distribution is also shown in Fig. 2(e). This distribution is strongly correlated 
with the photon energy distribution. Background was not subtracted in Fig. 2(e) because 
of the small number of signal and background events. Events mostly populate the low-mass 
range; however, the small data sample precludes any definite conclusions. It is clear that the 
observed (j)K~ mass distribution differs significantly from that in the three-body phase-space 
decay (also shown in Fig. 2(e)). To provide the correct branching fraction measurement, the 
MC data sample was adjusted to follow the measured (j)K~ invariant mass distribution. 

The Mbc distribution for the B^ (pK^'y decay channel is shown in Fig. 2(f). The figure 
is obtained in a similar manner to those for B^ (pK^'-f decay. In the fit of the Mbc 
distribution the peak position is fixed to the value 5.279 GeV/c^ and the background shape 
(defined by the slope parameter of the ARGUS function) is constrained to be flat. 

The signal yields, efficiencies, branching fractions and significances obtained for the B~ — >■ 
(j)K~'y and B^ (pK^'-f decay channels for the photon energy range 2.0 GeV < E* < 2.7 GeV 
are given in Table 1. The signal in the decay channel B~ (j)K~'y has a 5.5 a significance, 
whereas that for B^ — > (pK^'j is 3.3 a and a 90% confidence level upper limit, which is 
calculated assuming a Gaussian distribution of the statistical error and an additional one 
unit of the systematic error contribution, is also given. The statistical significance is defined 
as a/— 2 \n{Co/Cmax), where Cmax and Co are likelihood values at the best-fit signal yield 
and the signal yield fixed to zero. The branching fractions for the decays of the intermediate 
(p and states are taken from Ref. |17||] and are not included in the efficiencies quoted 
in Table 1. No event with more than one B candidate is found in the data; the double 
counting probability for the B~ (pK~'~f channel is estimated by MC to be (1.2 ± 0.3)% 
and is accounted for in the efficiency calculation. An equal production rate for the neutral 
and charged B mesons is assumed. 



TABLE I: The signal yields obtained from the Mbc fit, efficiencies, branching fractions, upper limit 
and significances for the B~ — > cpK~^ and — s- (pK^^ decay modes. 



Decay mode 


Yield 


Efficiency 


Branching fraction 


Significance 






(%) 


(10-6) 




B- (f)K--y 


21.6 ± 5.6 


12.9 


3.4 ± 0.9 ± 0.4 


5.5 




5.8 ± 3.0 


7.9 


4.6 ± 2.4 ± 0.6 


3.3 








< 8.3 (90% CL) 





In addition to the dominant continuum background, various BB background sources 
were studied. The only significant background found, with a contribution exceeding 1% of 
the signal level, arises from non-resonant B~ — > K~ K~ K^'^ decays. This contribution is 
estimated to be (4 ±4) % using events from the high mass (p meson sideband 1.05 GeV/c^ < 
Mk+k- < 1.25GeV/c^. The shape of this background as a function of M^+k- is taken 
from the MC simulation, assuming non-resonant decay i^'*-(1770) —>■ K~K~K^. The 
B~ — > (pK~'-f branching fraction is corrected for this contribution, and the uncertainty is 
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included in the systematic error. 

The major sources contributing to the systematic uncertainty of the branching fraction 
measurements are: the photon reconstruction efficiency and energy scale (4%), the recon- 
struction efficiency of charged tracks (1% per track), the charged kaon particle identification 
(1% per particle), the neutral kaon reconstruction efficiency (3%), uncertainties in the non- 
resonant background under the signal (4%), uncertainties in the angular distributions 
applied in the MC simulation of the hadronic states (5%), uncertainties in the MC simu- 
lation of tt" and rj rejection criteria (2%), uncertainty in the efficiency of the topological 
likelihood cut (2%), uncertainties of the background and B signal shape description in the 
fit (5%), and uncertainty in the determination of the number of BB pairs (< 1%). The sys- 
tematic uncertainties described above are added in quadrature to obtain the final systematic 
errors. 

In conclusion, the B ipK'-j decay modes were studied for the first time. A branching 
fraction for the B^ (f)K~'~f decay mode and an upper limit for the B^ (pK'^'y decay 
mode were obtained. The measured B^ —* (f)K~^ decay branching fraction is roughly an 
order of magnitude smaller than the B K*^ branching fractions. This factor can be 
attributed to the creation of the additional ss pair. Signal events are concentrated in the 
photon energy range 2.1 GeV < E* < 2.5 GeV, which is similar to the photon energy range 
for two-body radiative decays. 
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